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Abstract 

Scaling relations are one of the hallmarks of nodal superconductivity since they contain infor- 
mation characteristic for gapless order parameters. In this paper we derive the scaling relations for 
the thermodynamics and the thermal conductivity in the vortex state of the A and B phases of the 
skutterudite PrOs4Sbi2- Experimental verification of these scaling relations can provide further 
support for anisotropic gap functions which were previously considered for this material. 

PACS numbers: 74.25.Bt 
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1. Introduction 

Superconductivity in the filled skutterudite PrOs4Sbi2 was discovered in 2002 by Bauer 
et al [ll 0, [| , and has since generated ever- increasing attention. In particular, the presence 
of at least two distinct phases, the A and B phase, in an applied magnetic field is of great 
interest. Experimentally, it was observed that both phases have point nodes, and that the 
pairing channel appears to be a triplet with chiral symmetry breaking. ^,0,0]. However, 
the precise position of the A-B phase boundary is still controversial. For example, Measson 
et al found the A-B phase boundary to be almost parallel to H c2 (T) of the A phase. 
A possible explanation of this phase diagram was recently proposed in terms of the gap 



functions 



3 

( 

A B (k) = de^il-kt). (2) 



A A (k) = de^~(l- k$- k*- k% (1) 



Here e ±i(t>1 = (fc 2 ± ik 3 )/Jk% + fcf, e ±i(t>2 = (k 3 ± zJfei) / J k\ + k\, and e ±i(t>3 = {k x ± 



ik,2) / \J k\ + kl- The factor of 3/2 in the definition of A^(k) ensures proper normaliza- 
tion of the angular dependence of the order parameter. Furthermore, in Eq.(2) we choose 
the noda, d.tectton to be parallel to [001], because this p + h-wave order parameter symroetry 
is consistent with the magnetothermal conductivity data of Izawa et al pj. 

In 1997, Simon and Lee 0| introduced scaling relations for d-wave superconductors. More 



recently, following Volovik's approach 
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Kiibert and Hirschfeld 



llj obtained a scaling 



function for the quasiparticle density of states (DOS) in the vortex state of d-wave super- 



re thermodynamic 
From their very 



conductors. This expression for the DOS contains the scaling relations of t 
response tactions as well a s the thecal conduit, H^GH ^ 
general derivation it is clear that such scaling laws must apply to all nodal superconductors 
which have a comparable low-energy quasiparticle DOS G(E)~ \E\/A for \E\ < 0.3A in the 
absence of a magnetic field. If the above proposals (Eqs. (1) and (2)) for A(k) are correct, 
both phases of PrOs4Sbi2 would fall into this category. Experimentally, scaling laws for the 
specific heat have been verified experimentally in the cuprate superconductor YBCO [16], 
in the the ruthenate superconductor Sr2Ru04 with a magnetic field H || [001], and in 
the thermal conductivity of the heavy-fermion superconductor UPt3 |l£( . 

These measurements are consistent with the theory of scaling in nodal superconductors. 
Hence, scaling relations can be regarded as one of the hallmarks of nodal superconductivity. 



So far, however, scaling laws have only been studied in superconductors with line nodes, such 
as d-wave and f-wave order parameters. The object of this work is to extend these early 
analyses to superconductivity with point nodes by focusing on the skutterudite compound 
PrOs 4 Sb 12 . [rj. 

2. Quasiparticle Density of States 

Let us first consider the quasiparticle density of states in this compound, using the gap 
functions for the A and B phases given by Eqs. (1) and (2). In the a 
field, the low-energy quasiparticle DOS can then be approximated by 



)sence of a magnetic 



G A (E) = 1\E\/A, (3) 
Gb(E) = l\E\/A. (4) 

These equations are accurate in the low-energy regime E < 0.3A. Furthermore, in the 

vortex state the effect of the supercurrent can be introduced by letting E — ► E — v • q, where 

I I 

v • q denotes the Doppler shift. Following the derivation of Kiibert and Hirschfeld [ll( we 
obtain 
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i ^ 
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Ga(E, H) = I^^rinM^M (5) 

8=1 

G B (E,U) = ±e 3 g(E/e 3 ). (6) 



Here, the scaling function is given by 



g(s) = ^(i + ^) lfl >i (7) 



-VI - s 2 + —(1 + 2s 2 ) arcsin(s), s<l (8) 
4 4s _ 



and 



u , 

£i = -VeHsm6 u (9) 

sin^ = (1 -sin 2 9 cos 2 0) 1/2 (10) 

sinfl 2 = (1 - sin 2 sin 2 0) 1/2 and (11) 

sin # 3 = sin# (12) 

(9, <fi) are the angles indicating the direction of the applied magnetic field H. Note that 

in these equations the effect of impurity scattering is neglected. Therefore this result is 



valid only in the superclean limit, i.e., (TA) 1 / 2 < \E\,e < A, where T is the quasiparticle 

n 

scattering rate of the normal state. |2QJ In order to observe scaling behavior it thus appears 
necessary to have T < 0.01A. If such a sample is available the DOS obtained above should 
then be accessible by scanning tunneling microscope measurements. As seen from Eqs.(5) 
and (6), both Ga(E,H) and Gg(E, H) obey scaling laws. In particular the scaling law for 
Gb(E, H) is the same as in d-wave superconductors. 
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FIG. 1: The functions Ga(E, H) and Gb(E, H) along various directions of the applied magnetic 
field. 

In Fig. 1 G A (E,U) is shown for H || [111],H || [110] and H || [100] and G B (E,H) for 
H _L [001]. Along specific field directions we obtain 

Ga{E,H) = (yfe/ej , for H || [111], (13) 
= ^(g^E/e) + g{E/e)), for H || [110], (14) 
= ^9(E/e), for H || [110], (15) 



where e = \\feE. Note that G A {E, H) for H || [111] and H || [110] look very similar. Also 

due to the cubic symmetry of |A(k)| in the A phase (see the insert in Fig. 1) the cases 
H || [111], H || [—111] and H || [1 — 11] etc. are equivalent. 

In the B phase, the specific heat, the spin susceptibility, the superfluid density, and the 
nuclear spin lattice relaxation rate are then given by 

C s (T,H)/ ls T = ~f(T/e) (16) 

X,(T,H) = ~h(T/e) (17) 

p s ||(T,H)/p s (0,0) = i-^Lh(T/e) (18) 

TT\T,H)/T£ = (^)V(T/e), (19) 

where e = e 3 and p s \\ denotes the current parallel to the nodes (i.e. J || [001]). This 
expressions contain further scaling functions, 

^) 3 ^^)sech 2 g) (20) 

dsg(s)sec\i 2 (J^j and (21) 

dssg\s)^[^) (22) 

These expressions can be expanded in the low-temperature and high-temperature limits, 
with asymptotics given by 

f(T/e) = 1 + ^(T/e) 2 + . . . , for T/e « 1 (23) 

27C(3)T 3 , /0 , e T 1 

= -T^- + T- ln 2 ^ + • • • for - » 1 24 

4/T e 47r 1 e 



f{T/e) 


3 
2tt 2 


^) 


h(T/e) 


27\ 


POO 

'o 


J(T/e) 


2T, 


POO 
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h{T/e) = 1 + ^(T/e) 2 + . . . , for T/e « 1 (25) 

7rln2T tt e , / ,2T. 2 \ r T 

+ — -In l + (— ) 2 +..., for->l (26) 



2 e 32 T V r / 

tt 2 m 2 



J(T/e) = 1 + -) +...,forT/e«l (27) 



,7T, 



1 ,7r7\ 



? U ( T } J + ---' for 7 >>1 (28) 

These scaling functions are the same as in d-wave superconductors [14j and are shown in 
Fig. 2, where we introduced F — f — ^g^, K = h- and G = J - g(f) 2 . 
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FIG. 2: The scaling functions F(T/e), G(T/e), and K(T/e 



In analogy, for the A phase Eqs. (16), (17) and (18) are replaced by 21 1 



i 3 

C S (T,H)/ 75 T = ^/(T/e,) (29) 



3A . 

3 

3 A 



X.(T,H) = ^-J2eMT/e t ) (30) 



i=i 

3 



P4 (T, H)/p s (0, 0) = 1 - ^ ^ eMT/et) (31) 

3. Thermal Conductivity 

In order to determine the thermal conductivity it is necessary to include the effect of 
impurity scattering because unlike the thermodynamic response functions treated in the 
previous section, the scaling function for the thermal conductivity depends on the strength 
of the disorder, i.e. whether the impurity scattering is in the Born limit or the unitary limit 
jl^ . As we shall see below, the scaling functions for PrOs4Sbi2 turn out to be particularly 
simple if the heat current is parallel to a pair of point nodes. On the other hand, in the B 
phase the heat current has to be parallel to the nodal directions, in order to see an appreciable 



a. 



heat current. It appears that this condition is realized experimentally as reported in 
Otherwise the thermal conductivity would be much smaller since it vanishes like T(^) 2 as T 
approaches zero. In the A phase the heat current is always appreciable, although the thermal 
conductivity loses the cubic symmetry in the vortex state, unless H is directed along some 
symmetric direction (for example, H || [111], [—111], etc.). 



Following the derivation of Refs. |22j, |23( , the thermal conductivity is given by 

3n f°° 2 /z 2 h(u 1 U. 



AmT 2 



r^^/£^^\sech 2 (^/2T) (32) 

Jo \ r(^,H) / 



where 



If |£-vq| 2 -A 2 / 2 » 

2 V |(o3-v-q) 2 -A 2 / 21 ' (33) 



and 



f = ImV(a> - v ■ q) 2 - A 2 / 2 (34) 

Here < ... > denotes the averages over the Fermi surface and vortex lattice jisj . In the 
superclean limit ((rA) 1//2 < tu, |v • q| < A) Co is given by 

u> = u + iT{— l^~ v "^l \ (35) 

\ v /(^-vq) 2 -A 2 / 2 / 

~ cj + 2rG(cj,H) (36) 

in the Born limit. And in the unitary limit we find 

u = uj + iTG- 1 (uj,B.) (37) 

where G(u, H) for the A and B phase have been defined in Eqs. (5) and (6). 

Let us first consider the Born limit. Substituting Eq.(35) into Eq.(32) we obtain 

fi a = / rfcucu 2 ^ — j — sech (u;/2T) 38 

8mT 2 r J < \lj - v • q| > 

where < \uj — v • q| >' is the same as < \u — v • q| > but with contributions from nodes 
at (001) and (00 -1) only. Then in the B phase < \lu — v • q| >'=< \lu — v ■ q| > and 
thus k^ z = |ft n , where n n = is the thermal conductivity in the normal state. The 
thermal conductivity in the Born limit is independent of H. On the other hand in the A 
phase < | a; — v • q| >V< |w — v • q| >, leading to 

ni/n n = i for T > e, (39) 
= ^e 3 /(ei + £2 + e 3 ), for T < e. (40) 

The e^'s were defined in Eq.(9). 
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The scaling law is of greater interest in the unitary limit. First let us consider the B 
phase where < \uo — v • q| >'=< \uo — v ■ q| >. Substituting Eq. (37) into Eq. (32) one 
obtains 

o 3ir 2 n 



512mr(TA) 2 

3ne 5 



duu 2 < \uj - v • q| > 2 sech 2 (w/2T) 



32mr(TA) 2 
and the scaling function is defined as 

k£(T,H) _ 120, T. 6 



poo 

\ dss 2 g 2 (s)sech 2 (es/2T) 
Jo 



F(T '^ " = r dssV(s)sech 2 (es/2T) 

K zz\ 1 ^ V ) " e JO 



where e = e%. 
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FIG. 3: The scaling function F(T/e) for the unitary limit (U), the Born limit (B) and the case 
without inversion symmetry (I) are shown as a function of T/e. 



where e = £3. This scaling function is shown in Fig. 3, with asymptotics given by 



and 



2l7r 4 / T, 



640 V A 
3tt 2 - 



1 + ?(^ 2 + 



8 V A 



(ir 1 



7tt 2 ,T 



15 v e 



-) 2 + 



, for e < T 
. ) , fore>T 



F B (T/e) = 



l + ^)V..,for e «T 
80 



7vr 4V A 



1 + 7F<7> + 



, fore > T 



(44) 
(45) 

(46) 
(47) 
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This scaling function F(T/e) is the same in other nodal superconductors such as those 
with d-wave symmetry. For exam ple, F B (T/e) describes very well the scaling behavior 
recently observed by Suderow et al 18[ in UPt3. 

In the A phase the scaling function is somewhat more complicated. We find 



. ne 3 



72mr(AT) 2 

In particular for H || [111] and [100] Eq.(48) reduces to 



/oo / 3 \ 

dww 2 g(u/e 3 ) f ^e^(o;/e 3 )J sech 2 (cu/2T) (48) 



«4 = 24mr( AT) 2 J Q dss2 9 2 (s)sech 2 (e 3 s/2T) forH || [111], and (49) 

POO 

/ dss 2 # 2 (s)sech 2 (e 3 s/2T) forH || [100] (50) 
Jo 



36mr(AT)' 



where e 3 = 2 e y g and 2 e for H || [111] and H || [100] respectively. Therefore in 
these two cases we will have the same scaling function as F B (T/e). 
4. Concluding Remarks 

We conclude that the scaling behavior of the universal heat conduction and the thermal 

n 

conductivity can be regarded as a hallmark of nodal superconductivity [20]. Moreover, the 
scaling function F B (T/e) describes the thermal conductivity data measured in UPt3 by 



Suderow et al 



18| very well. In this paper, we have found that the thermal conductivity in 



both the A and B phases of PrOs4Sb 12 exhibits a number of characteristic scaling relations. 
The directional dependence of these scaling relations on H and q is expected to further 
confirm the nodal structure of A(k) proposed in 

In the course of the present study we have also observed that the scaling behavior of the 
thermal conductivity in CePt 3 Si found by Izawa et al (24! 25] is very unusual. Their data 
appears to be more consistent with the case where the inversion symmetry of the impurity 
scattering is broken. Clearly, further study of scaling laws in nodal superconductors will 
open a new point of view on the whole subject. 
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